Non-conformal asymptotic behavior of the time-dependent field-field correlators of ID 

anyons 
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The exact large time and distance behavior of the field-field correlators has been computed for one- 
dimensional impenetrable anyons at finite temperatures. The result reproduces known asymptotics 
for impenetrable bosons and free fermions in the appropriate limits of the statistics parameter. 
The obtained asymptotic behavior of the correlators is dominated by the singularity in the spectral 
density of the quasiparticle states at the bottom of the band, and differs from the predictions of the 
conformal field theory. One can argue, however, that the anyonic response to the low-energy probes 
is still determined by the conformal terms in the asymptotic expansion. 
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When confined to move in one dimension, hard-core 
impenetrable particles cannot be exchanged directly. 
This mean that the symmetry of the wavefunction with 
respect to exchanges of coordinates can be defined ar- 
bitrarily, e.g., independently of the dynamic particle- 
particle interaction (see, e.g., [J, 0). The wavefunc- 
tion describes then the exchange statistics with, in gen- 
eral, fractional statistics parameter k € [0, 1] interpolat- 
ing between Bose (k = 0) and Fermi (k = 1) particles. 
Quantitatively, this "anyonic" exchange statistics iniplies 
that the multi-particle wavefunctions behave as 0, 0, [1] : 



upon exchange of the two nearest-neighbor coordinates, 
where e{x) = x/\x\, e(0) = 0. While the wavefunctions ip 
depend on the statistics parameter k, the absence of real 
exchanges means that thermodynamic properties of the 
hard-core particles are independent of statistics. In par- 
ticular, in the simplest case of only hard-core interaction 
of zero radius, which is considered in this work, the ther- 
modynamics coincides @, 0] with that of free fermions. 

The exchange statistics manifests itself, however, in 
the correlation functions of the field operators, which de- 
scribe processes with changing number of particles. In 
the general case of ID anyons, only "static" same-time 
correlators of the field operators have been calculated 
[1, 1^ [13, U, 121 so far. The purpose of this work is to 
present the exact calculation of the large time and dis- 
tance asymptotics of the time-dependent field-field corre- 
lators of impenetrable anyons. The time-dependent cor- 
relators are more relevant experimentally, as they deter- 
mine, e.g., the tunneling density of states of the system. 
For anyons, tunneling characteristics can be measured, 
e.g., in FQHE structures (see, e.g., dill). The new 
qualitative feature of the time-dependent anyonic corre- 
lators obtained in this work is that, in contrast to previ- 
ously found lj,ll5| asymptotics for impenetrable bosons, 
the correlators for the general statistics parameter k do 
not exhibit conformal behavior in the low-temperature 
limit. This result challenges the accepted notion that 



the large-time asymptotic behavior of the correlators of 
the ID systems is determined by the low-energy excita- 
tions. Our exact calculation shows that the large-time 
asymptotic behavior can actually be dominated by the 
singularity of the density of the quasiparticle states at 
the bottom of the energy band which is characteristic 
for one dimension. This conclusion adds another limita- 
tion on the applicability of the conformal field theory to 
the large space and time asymptotic description of the 
ID particle systems. Previously known limitations were 
related to the effects of the strong interaction on parti- 
cles with spin, which leads to vanishing of the exchange 
energy scale, making the spin dynamics non-conformal 
IgL 17 , [H, |3| , and the effects of finite nonlinearity of 
the quasiparticle spectrum at the Fermi energy ji^l . 

In details, the model considered in this work describes 
one-dimensional hard-core anyons interacting via a S- 
function potential, which simulates the hard-core re- 
pulsion, and are free otherwise. The second quantized 
Hamiltonian of such anyons is given by 

-/i*^(a;)«'(a;)] , (1) 

in the c — *■ cxd limit. Here h is the chemical potential, 
and the anyonic fields satisfy the following commutation 
relations 

^(a;i)*^(a;2) = e-'"«^("i-"^)^t(2;2)^(a;i) -f ,5(xi - xa), 
^'(a:i)^'(a;2) = e'^''^(^i^^^)^'(a;2)*(2:i) . (2) 

The focus of this work is on the calculation of the large 
time and distance asymptotic behavior of the field-field 
correlation function at finite temperature defined as 

As a function of the statistics parameter k, this correla- 
tor should interpolate between the case of impenetrable 
bosons (k = 0) and free fermions {k ~ 1). 
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To state first our main result, we introduce the rescaled 
variables x = (xi - X2)Vt/2 > 0, t = {t2 - ti)T/2 > 
0, (3 — h/T. In terms of tliese variables, the field corre- 
lator can be written as 



bosons and fermions. In the bosonic limit k 



(*(a;2, i2)*^(a;i, ii)) = VTg{x, t, /3, k) 



(3) 



where the function g is defined below. Due to the fact 
that g{x,t, 13, k) — g{—x,t,P,—K) and g{x,t, fi, k) = 
g*(x, —t, P, —k) it is sufficient to consider only the case 
a; > 0, t > 0. One needs to distinguish the space-like, 
x/2t > y/p, and time- like, x/2t < ^/]3, regions. While 
the leading term in the asymptotics is the same in both 
regions, the next-to-leading term behaves in them differ- 
ently. The large time and distance asymptotic form of 
the function g, for x,t oo, with x/t = const, is: 



g{x,t,l3,K) 



^ly^ /2^C(x,t,P,K)+ixI(l3,K} 



(4) 



where 

AT ^ [_ai/2^,(^_ 2/5)1/2]/^^ a = 13+^(3^ + n^K^, 

and the upper (lower) sign corresponds to the space- 
like (time-like) regions. Also, — —x/2t, v = 
— i In |(^(As, /3, k)|, the constants Cq, C\ are some undeter- 
mined amplitudes, /(/3, k) = Q ln(p(A2,/3,K) dA/TT, 
and 



C[x, t, (3, k) — 



1 



2U|ln|(^(A^/3,K)| dX, (5) 



with 



^{X\(3, k) ^ {e^"-f' - e^^^)/{e^'-^ + 1) . (6) 



Even though it seems that the second term (with ampli- 
tude ci) in the expansion ([4]) is superfluous, because it is 
in general exponentially decreasing and smaller than the 
error term, one can see that it gives the dominant con- 
tribtion in the bosonic limit, thus justifying its presence. 
Another remark is that the precise analytic expression for 
the second term in the expansion (|4]) depends on the an- 
alytical structure of the solution of the Riemann-Hilbcrt 
problem. While there is no additional difficulties in an- 
alyzing all regimes, Eq. ^ is written out specifically in 
the situation when this term is determined by the pole 



at , which is the case when 



|5RV/3 + «7rK - x/2t\ > 3/9 + 



(7) 



The positive branch of the square root should be taken 
in this relation. This condition is always satisfied, in par- 
ticular, in the more interesting low-temperature regime, 
(3 » 1. 

We begin the analysis of Eq. ^ by demonstrat- 
ing that it reproduces the known asymptotics for 



we have ip{X'^,P,0) = (e^ "'^ 
-(l/7r)ln|(p(A2,/?,0)|, and 

C{x,t,l3,0) = - / |x- 2U|ln 



sA2-/3 + 1 



1), 1^ 



dX. 



Also, /(/3,0) = -2y^ for /3 > 0, and /(/3,0) = for 
/3 < 0. In the case of negative chemical potential, /3 < 0, 
condition ([7]) implies that the second term in the expan- 
sion ([4]) is exponentially small and the asymptotic be- 
havior is determined by the first term, which reduces for 
K = to: 

g{x,t,P) ~ cot(''-')'/2gC(.,t,/3,0)+2.t(A^+/3)^ 

For the positive chemical potential, however, A J = — vT^, 
and the leading contribution in both space-like and time- 
like regions is given by the second term in the expansion 
O, thus giving for f3 ~ h/T > 0: 

gix,t,(3)^cy/^e^^-''-'^''>\ 

Both of these expressions agree with the previous results 
14, 1^ for the impenetrable bosons. 

In the case of free fermions, the equivalent g function 
in the rescaled variables is 



g{x,t,(3) = 



„2it/3 

"2^ 



+00 gA^-/3 



-2itA^-2ia;A 



dX. 



In the large t and x limit with t/x = const, the steepest 
descent method gives the following asymptotic behavior 
under the condition 

g{x,t,P)~cot-'/^e^''^^+^^^ +cie^^^'''-"^f\ 

Here Xj are given by the same expression as Aq^ with 
K = 1, and as before, the upper (lower) sign corresponds 
to the space-like (time-like) region. Comparing this ex- 
pression to Eq. ([4]), which is simplified for k — > 1 by 
the relations C{x,t,(3, 1) — I{[3, 1) = v ^ Q, we see that 
Eq. (ID) indeed reproduces the correct behavior of the 
asymptotic correlators for free fermions. 

At positive chemical potential and low temperatures 
the system is expected to exhibit conformal behavior. 
One of the simplest ways to calculate the asymptotic cor- 
relators in this reg ime is to use the standard bosoniza- 
tion approach 2l|, in which the field operators are ex- 
pressed through two bosonic fields: the integral 9{x) of 
the long- wave part of the density fluctuations p(x), i.e., 
d9{x)/dx = p{x)/tt, and the conjugate fleld (j){x) defined 
by the commutation relation [(f){x),9{x')] = iTre{x—x')/2. 
The fact that the operator e*"^ reduces density by the 
amount that corresponds to one particle, and 6 produces 
the appropriate phase changes of the wavefunction across 
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each particle, implies then that the anyonic operators ^ 
can be written as [231 



*(x,t) 



^-i{K+2m)[kpx+e{x,t)] ^i. 



{x,t) 



(8) 



The sign ~ in Eq. ([5]) is the reminder that the relative 
amplitudes of different components are not defined in 
the sum. The commutation relations between and 
shows directly that the individual terms in Q satisfy the 
appropriate exchange relations ^ away from xi — X2- 
The standard average over the equilibrium fluctuations of 
9 and (f> gives then the large time and distance behavior 
of the correlator of the anyonic field operators ([8]) 0, H^l • 
The two leading asymptotic terms are: 



(^'(x2,i2)«'Ua;i,ti)) 



bo e*'''=^^i^ 



l(l-K)V4f,, ,1(1+«)V4 



+ 



b-i e^(''-2)fcF2;i 



[zi_](3-«)V4[y^](l— )V4 



(9) 



where u± = sinh[7rT(i2i — iO ± xi2/vf)], the constants 
b are some undetermined amplitudes of different compo- 
nents of the correlator, and i2i = t2 — ii, 2:12 — xi — X2- 
The Fermi vector and the Fermi velocity are, respectively, 
kp = Vh and vp = 2\fh. 

In Equation ^ , we kept the two leading terms to facil- 
itate comparison to the case of the static correlators [13] , 
for which both terms are important in the vicinity of the 
Fermi point k — \. For the dynamic correlators consid- 
ered in this work, it is assumed that t2i 7^ 0. In this 
case, only the first term is relevant for our discussion. In 
the T ^ limit, the criterion used above to separate the 
space-time and time-like regions reduces to the compar- 
ison between ^21 and xyilvp. Keeping only the leading 
term in Eq. ((H) we see that this equation reduces to 



(*(x2,t2)4'^(a^i,ti)) 



in the space-like region, and to 

(vl/(x2,i2)*^(xi,ii)) ~ e*fc-''^"e'^"^^e-""^^(^+'^), 

in the time-like region. In the same low-temperature 
limit, one can also simplify our main expression In- 
deed, for T ^ one has: ixl{(3, k) ~ 2ix\/P{K — 1), and 
X^-^^Ti^- Also, 



C{x,t,fi,K) = -7r(l - nf 



xl{2^), x/2t>^, 
t, x/2t<^/]3. 



Using these equations one can see directly that the ex- 
ponential parts of the next-to-leading term in Eq. (|4]) 
combine to give exactly the same result as predicted by 
CFT. This means that although the leading asymptotic 
term in the correlator ^ is not conformal, the next-to- 
leading term is. The exponential factors in the leading 



term of the correlator (|4]) imply that in the low temper- 
ature limit, /3 ^ 1, this term can manifest itself in the 
response of the anyonic liquid to external perturbation 
only in the case of large energy u; of the perturbation 
Lo ep = h. Therefore, at small energies, u <C cf, the 
response is determined by the second asymptotic term, 
and coincides with the prediction of the CFT. 

To see this more explicitly, one can consider the sim- 
plest example of weak single-point tunneling into the any- 
onic liquid. As usual, the tunneling density of states is 
given in this case by the Fourier transform of the field 
correlator with x — Q. Equation (jl]) for the anyonic cor- 
relator reduces for a; = (as before, /? ^ 1) to: 



cie 



(10) 



Taking Fourier transform of this correlator, one can see 
that the contribution of the first term to the tunneling 
density of states Alto) is A{uj) cx (ft, — uj)~'^/'^. Although 
the large-time asymptotics ([TUl) is not sufficient to es- 
tablish all features of Aiuj), this estimate illustrates the 
fact that this term in the asymptotic correlator is pro- 
duced by the singularity of the ID density of states at 
the bottom of the band. This implies that at small en- 
ergies this contribution to tunneling will be suppressed, 
and the tunneling density of states will exhibit the power- 
law nonlinearities with the power (l-f k^)/2 produced by 
the leading term in the CFT correlator 

As the last part of the discussion of our main result 
(|4]), we provide a simple heuristic interpretation of the 
exponential factors in the leading term in this equation. 
The qualitative statement proven by this interpretation is 
that in the T 7^ regime considered in this work, the ex- 
ponential decay of the correlator is caused mainly by the 
thermal fluctuations of the number of particles which cre- 
ate the fluctuations of the accumulated statistical phase. 
The average taken over these phase fluctuations leads to 
the exponential suppression of the correlator. To see this, 
one starts with the standard Wigner-Jordan transforma- 
tion expressing the anyonic fields in terms of (in our 
case free) fermions ^: 

*(x,t) =e*'^(i-'")"(="'*)e(a:,t), n{x,t) = ( dx'p{x',t), 



where p = is the operator of the particle density. This 
transformation results in the following expression for the 
anyonic field-field correlator: 

{^ix2,h)^\xi,h)) = (e(x2,t2)e*"("-i)"^e^(a;i,ii)), 

where rid = n{xi,ti) — n{x2,t2)- For large distances 
xi — X2 and time differences t2 — ti, the fluctuations 
of Ud can be calculated from the following quasiclassical 
considerations. Contribution to these fluctuation from 
particles with different momenta k can be treated in- 
dependently from each other. For given k, the particle 
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fluctuations at the initial moment ti are governed by the 
Fermi distribution = l/[e^'^ -'»)/^ -|- i]^ which gives 
the probabihty density 'd{k)dk/2Tr for momentum-fc par- 
ticles to be present in any space interval. The resulting 
random distribution of particles moves with time with ve- 
locity (in our notations) 2k. The dilTerence nd can then 
take values 1,0, or -1, producing the phase factors in the 
correlator — e*^**, 1,— e^*'^''. Analyzing the shifts of the 
distributions with time, and averaging over the proba- 
bilities to have or not to have a particle in the initial 
distribution, one sees that 



B++ = e+(A)(l + VT)"^e+(A) dX. The integral op- 



^C(x,t,f3,K)+ixI{p,K) 



(11) 



demonstrating that this part of the asymptotic suppres- 
sion of the anyonic correlator is indeed caused by thermal 
fluctuations of the statistical phase. 

Qualitatively, the non-conformal contribution to the 
asymptotic behavior of the anyonic field-field correlator 
found in this work is analogous to the power-law correc- 
tions to the exponential decay of the metastable states 
in general quantum mechanical problems (see, e.g., 
and references therein). Similarly to these general correc- 
tions [23| , which also result from the bounded- from-below 
nature of the energy spectrum of a quantum system, 
our non-conformal power-law terms in the asymptotics 
come with the energy shift to the bottom of the band. 
This shift ensures that despite the presence of the non- 
conformal terms, the anyonic response at energies close 
to the Fermi energy, agrees, as demonstrated explicitly 
by our results, with the predictions of the conformal field 
theory. 

Finally, we present a brief outline of the derivation of 
the asymptotic behavior of the correlator ^ . The start- 
ing point of our calculations is the determinant repre- 
sentation of the space-, time-, and temperature depen- 
dent field-field correlator for anyons obtained in In 
rescaled variables, the correlator is given by Eq. with 



]{x,t,l3, k) 



27r ^ 



.det{l + VT)- 



(12) 



In Equation p2)) . Vt is an integral operator acting on 
the entire real axis with the kernel given by 



VT{x,^l) 



e+{X)e_{n) - e_(A)e+(/j) 
X — fj. 



where 



e-(A) 



cos(7rK/2) 



, e+(A) = e_(A)ii;(A) 



E{X) = P.V. 



-2i0(Ai) 

dfi -I- 7rtan(7rK/2)e^2"^(^' 

^ - A 



and (f>{X) = tX^ + a;A. The auxiliary potential &_)__|_ is de- 



erator appearing in Ecl (1121) is of the special kind called 
"integrable" operators[15, 25]. This type of operators ap- 
pears frequently in the investigations of correlation func- 
tions of integrable models and distribution of eigenvalues 
of random matrices 25| . The specific factorization prop- 
erties of the kernel Vt allow one to obtain differential 
equations for the correlators. These differential equa- 
tions do not depend on the statistics parameter k and are 
the same as the ones obtained for impenetrable bosons. 
The statistics parameter enters in the initial conditions 

[l2| . The large time and distance asymptotic behavior 
of the correlation functions is obtained from the analy- 
sis of a matrix Riemann-Hilbert problem associated with 
these differential equations using the "Manakov ansatz" 
11, [11 EE [13 . The details of the computations will be 
presented elsewhere. 

The authors were supported by the NSF grants PHY- 
0653342, DMS-0503712 and DMR-0325551 . 



fined as b++ = B++ - G with G = / 



+ 00 



-2i0(A) 
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